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Time-dependent density functional theory (TDDFT) is presently enjoying enormous popularity 
in quantum chemistry, as a useful tool for extracting electronic excited state energies. This article 
discusses how TDDFT is much broader in scope, and yields predictions for many more properties. 
We discuss some of the challenges involved in making accurate predictions for these properties. 



Kohn-Sham density functional theory 0, 0, Q is the 
method of choice to calculate ground-state properties of 
large molecules, because it replaces the interacting many- 
electron problem with an effective single-particle problem 
that can be solved much faster. Time-dependent density 
functional theory (TDDFT) applies the same philosophy 
to time-dependent problems. We replace the complicated 
many-body time-dependent Schrodinger equation by a 
set of time-dependent single-particle equations whose or- 
bitals yield the same time-dependent density n(r, t) . We 
can do this because the Runge-Gross theorem 4] proves 
that, for a given initial wavefunction, particle statistics 
and interaction, a given time-dependent density n(r,t) 
can arise from at most one time-dependent external po- 
tential «ext(r,t). We define time-dependent Kohn-Sham 
(TDKS) equations that describe TV non-interacting elec- 
trons that evolve in V s (r, t), but produce the same n(r, t) 
as that of the interacting system of interest. Develop- 
ment and applications of TDDFT have enjoyed exponen- 
tial growth in the last few years |f| 0, Q- an d we hope 
this merry trend will continue. 

The scheme yields predictions for a huge variety of phe- 
nomena, that can largely be classified into three groups: 
(i) the non-perturbative regime, with systems in laser 
fields so intense that perturbation theory fails, (ii) the lin- 
ear (and higher-order) regime, which yields the usual op- 
tical response and electronic transitions, and (iii) back to 
the ground-state, where the fluctuation-dissipation theo- 
rem produces ground-state approximations from TDDFT 
treatments of excitations. 

In the first, non-perturbative regime, we have sys- 
tems in intense laser fields with electric field strengths 
that are comparable to or even exceed the attractive 
Coulomb field of the nuclei 5] . The time-dependent 
field cannot be treated perturbatively, and even solving 
the time-dependent Schrodinger equation for the evolu- 
tion of two interacting electrons is barely feasible with 
present-day computer technology 0- For more electrons 
in a time-dependent field, wavefunction methods are pro- 
hibitive, and in the regime of (not too high) laser intensi- 
ties, where the electron-electron interaction is still of im- 
portance TDDFT is essentially the only practical scheme 



available. With the recent advent of atto-second laser 
pulses, the electronic time-scale has become accessible. 
Theoretical tools to analyze the dynamics of excitation 
processes on the attosecond time scale will become more 
and more important. An example of such a tool is the 
time-dependent electron localization function (TDELF) 
|ToL fTl| - This quantity allows the time-resolved obser- 
vation of the formation, modulation, and breaking 
of chemical bonds, thus providing a visual understand- 
ing of the dynamics of excited electrons (for an example 
see FigQ]and Ref. 01 )■ The natural way of calculating 
the TDELF is from the TDKS orbitals. Recent appli- 
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FIG. 1: Snapshots of the time-dependent ELF for the excita- 
tion of acetylene by a 17.5 eV laser pulse |T^|. The pulse had 
a total length of 7 fs, an intensity of 1.2 x 10 14 Wcm -2 , and 
was polarized along the molecular axis. Ionization and the 
transition from the bonding n state to the anti-bonding ir* 
state are clearly visible. 



cations in the beyond perturbative regime range from 
above-threshold ionization of metal clusters 0] to coher- 
ent control of quantum wells|l4| to multiharmonic gen- 
eration in benzene|l5l|. 

A much larger group of applications in chemistry is the 
linear response to a spatially uniform electric field, i.e., 
the optical response in the dipole approximation | la 
Kl\ . Formal analysis of this situation shows that TDDFT 
yields predictions for electronic excitations, both their 
position (transition frequency) and intensity (oscillator 
strength). These are corrections to transitions between 
occupied and unoccupied levels of the ground-state KS 
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potential, thus providing a simple interpretation of those 
levels [3. In the area of calculating electronic excita- 
tions, TDDFT is rapidly becoming a standard tool, com- 
plimentary to existing wavefunction techniques^! ■ Just 
as in the ground-state case, it has the advantage in com- 
putational speed, allowing study of larger systems than 
with traditional methods, and the usual disadvantage (or 
excitement) of being unsystematic and artful. A final 
application is to write the ground-state XC energy 
in terms of the frequency-dependent response function, 
and so linear response TDDFT yields approximate treat- 
ments of the ground-state problem |2fJ. I2H l22l |23| . 

A random walk through some of 2004's papers using 
TDDFT gives some feeling for the breadth of applica- 
tions. Most are in the linear response regime. In in- 
organic chemistry, the optical res pon se of many transi- 
tion metal complexes (H HHH |H |H H 
calculated, and even some X-ray absorption 
ganic chemistry, the response of thiouracil 
tetrazine0], an d annulated porphyrins 01 were inves- 
tigated. In photobiology, potential energy curves for 
the trans-cis pho to-isomerization of protonated Schiff 
base of retinal[3^| have been calculated. For these and 
other systems, there is great interest in charge-transfer 
excitations jH HH HI E| HI H| , but (as we later dis- 
cuss) intermolecular charge transfer is a demanding prob- 
lem for TDDFT. Another major area of application is 
clusters, large and small, covalent and metallic, and ev- 
erything inbetween[4ll Hi Q EE! E3 , including 
Mct-Cars^. Several studies include solvation, for ex- 
ample, the behavior of metal ions in explicit water[4flj|. 
TDDFT in linear response can also be used to calculate 
both electronic and magnetic circular dichroism|50l l5l| . 
and has been applied to helical aromatics[52j. and to 
artemisinin complexes in soluti onl53l. T here have also 
been applications in materials |54l Ififij and quantum 
dotsjE3 but, as discussed below, the optical response of 
solids requires some non-local approximations 57J. Be- 
yond the linear regime, there is also grow in g^ interest in 
second- and third-order response |58l l59l l&oL l6l| in all 
these fields. 

A wonderful aspect of TDDFT is that a single ap- 
proximation to the time-dependent XC potential implies 
predictions for all these quantities. This is analogous 
to the ground-state case, where a single approximation 
to E xc can be applied to all ground-state properties 
of all electronic systems, such as dissociation energies, 
bond lengths and angles, vibrational frequencies, etc., of 
atoms, molecules, clusters, and solids. The starting point 
of most TDDFT approximations is the adiabatic local 
density approximation (ALDA), which approximates the 
XC potential at point r and time t by that of a ground- 
state uniform electron gas of density n(rt). This is clearly 
accurate when the density varies sufficiently slowly in 
time and space, but works surprisingly well beyond that 
limit for many systems and properties, just as LDA does 



for most ground-state problems. 

We make an important distinction here between the 
mature subject of ground-state DFT, and the developing 
one of TDDFT. In the former, there is general consen- 
sus on which properties are captured by which function- 
als, and the aim is toward higher accuracy^- One ex- 
pects chemical bonds to form in modern KS DFT calcula- 
tions, and one hopes to use better functionals to produce 
better accuracy and reliability • But time-dependent 
quantum mechanics probes a far more diverse range of 
electronic phenomena, and in TDDFT, we are still ex- 
ploring even which properties are captured at all by the 
presently available approximate functionals. Quantita- 
tive accuracy is less of an issue as yet. Most data on the 
performance of TDDFT are for systems driven by some 
external field. Practically nothing is known about how 
TDDFT performs in the description of relaxation pro- 
cesses, i.e., on the time evolution of large systems starting 
from a non-equilibrium initial condition j64|. A closely re- 
lated question is the description of density fluctuations 
within TDDFTjsU. 

At this point, we introduce a few equations, to make 
the discussion more precise. We use atomic units 
throughout, and suppress spin indices. For brevity, we 
drop commas between arguments wherever the meaning 
is clear. We write the TDKS equations as 



. d4> 3 {rt) 
1 dt 



■w s [n](rt)Uj(rt) 



(1) 



whose density n(rt) = J2j=i l<Aj( r OI is precisely that 
of the real system. We define the exchange-correlation 
potential via 

r n(r't) 

v s (rt) = zw(ri) + / dV + v xc (rt). (2) 

J |r-r'| 

The exchange-correlation potential, v xc (rt) is in general 
a functional of the entire history of the density, n(rt), 
the initial interacting wavefunction \&(0), and the ini- 
tial Kohn-Sham wavefunction, $(0). This functional is 
a very complex one, much more so than the ground- 
state case. Knowledge of it implies solution of all time- 
dependent Coulomb interacting problems. If we always 
begin in a non-degenerate ground state|66|. the initial- 
state dependence can be subsumed by the Hohenberg- 
Kohn theorem^]], and the only unknown part of v s (rt), 
the exchange-correlation (XC) potential, is a functional 
of n(rt) alone. 

In the special case of the response of the ground state 
to a weak external field, the system's response is charac- 
terized by the non-local susceptibility 



dt' / d 3 r' X [no}(r,r';t~t') Sv ext (r't'). (3) 



Sn(rt) 



X is a functional of the ground-state density, no(r). 
The central equation of TDDFT linear response [lg is 
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, Dyson-like equation for the true x °f the system: 

X(rr'cj) = x s (rr'w)+ [ d 3 n [ d 3 r 2 Xs(rriw) 
1 



ri -r 2 



+ / xc (r ir2 cj) x(r 2 r'w), (4) 



Here Xs is the Kohn-Sham response function, constructed 
from KS energies and orbitals: 



$„(r) $*(r') 



where g is a double index, representing a transition from 
occupied KS orbital i to unoccupied KS orbital a, uj q — 
e a - a, and $,(r) = 4>* (r)0„(r). Thus Xs is purely a 
product of the ground-state KS calculation. On the other 
hand, the XC kernel is defined as 

/ xc [n ](rr',i-t') = Sv xa (rt)/dn(r't')\ no . (6) 

This is a much simpler quantity than v xc [n] (ri), since the 
functional is only evaluated at the ground-state density. 
It is non-local in both space and time. The non-locality 
in time manifests itself as a frequency dependence in the 
Fourier transform, / xc (rr'w). 

Next, we introduce Casida's equations [12], in which 
the poles of x are found as the solution to an eigenvalue 
problem: 



9' 



^Vl qql (uj) a q i = lo 2 a q , 



(7) 



where 



Q qql (uj) = u; q 6 qq > + 2^Ju q ~u/ q (q\f lixc (uj)\q'}. (8) 

and (<7|/ H xc(k')|<7') is the matrix element of the (Hartree)- 
XC kernel in the set of functions $ 9 (r). Eigenvalues yield 
the square of transition frequencies, while eigenvectors 
yield oscillator strengths. Ignoring off-diagonal matrix 
elements can yield much insight into the nature of the 
TDDFT corrections to the KS transitions [l^. 

Lastly, we mention how TDDFT produces sophis- 
ticated approximations to the ground-state exchange- 
correlation energy. The adiabatic connection fluctuation- 
dissipation formula is: 

2W«o] = -~ [dX J >r / 'dV ^ 

^ {x«N(rr'o;) + n (r)<5( 3 '(r - r')](9) 

where the coupling-constant A is defined to multiply 
the electron-electron repulsion in the Hamiltonian, but 
the external potential is adjusted to keep the density 
fixed |r3^. |68| . So any model for / xc , even setting it to 
zero (called the Random Phase Approximation) , yields a 



sophisticated approximation to E xc , by solving Eq. Q 
for x (at each A) and inserting in Eq. (JSJ. 

All the above equations are formally exact. In 
any practical DFT calculation, approximations must be 
made. The most common approximation in TDDFT is 
the adiabatic approximation, in which 



<f[n](ri) = < c [n Q ](r)|„ o(r)=n(rt) , 



(10) 



i.e., the XC potential at any time depends only on 
the density at that time, not on its entire history. 
This becomes exact for slow variations in time. Most 
applications, however, are not in this slowly varying 
regime. Nevertheless, results obtained within the adi- 
abatic approximation are, in most cases, rather accu- 
rate. Any ground-state approximation (LDA, GGA, hy- 
brid) automatically provides an adiabatic approximation 
(e.g., ALDA) in TDDFT. Moreover, the XC kernel is 
frequency-independent in the adiabatic approximation, 
taking its u> — > value. 

As mentioned above, TDDFT is proving very useful 
in predicting optical response properties of molecules. 
The Casida equations have been encoded in most stan- 
dard quantum chemical packages, and efficient algo- 
rithms developed to extract the lowest-lying excitations. 
A small survey is given by Furche and Ahlrichs 19] . Typ- 
ical chemical calculations are done with the B3LYP |69j 
functional, and typical results are transition frequencies 
within 0.4 eV of experiment, and structural properties 
of excited states are almost as good as those of ground- 
state calculations (bond lengths to within 1%, dipole mo- 
ments to within 5%, vibrational frequencies to within 
5%). Most importantly, this level of accuracy appears 
sufficient in most cases to qualitatively identify the na- 
ture of the most intense transitions, often debunking 
cruder models that have been used for interpretation for 
decades. This is proving especially useful for the pho- 
tochemistry of biological molecules [jjj. An alternative 
implementation, often favored by physicists, is to propa- 
gate the TDKS equations in real time, having given the 
system an initial weak perturbation. Such calculations 
either use a real-space- grid fFH Fizj or plane waves [7^. 

This article is not about the (admittedly) gratifying 
successes of TDDFT calculations, which are dicussed in 
recent reviewsj^, @ and the recent literature. We be- 
gin from there, and explore a much wider arena. To do 
this, in Fig. [3 we have drawn a cartoon (literally, a stick 
figure) to represent the information in a typical calcula- 
tion. Each line represents a transition, with its position 
denoting the transition frequency and its height propor- 
tional to the oscillator strength. For the He atom, the 
2s— >2p singlet transition is at 21.2 eV, while for the N 2 
molecule, the 3 n 7 transition is at 7.4 eV. Many applica- 
tions of TDDFT report only the positions of the few low- 
est optically-allowed transitions, while some report also 
their oscillator strengths. 
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FIG. 2: Cartoon of the exact optical absorption spectrum of 
an atom or molecule, with discrete transitions represented by 
straight lines (see text). 
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FIG. 3: Same as Fig. [5] but showing higher frequencies, 
including the infinite Rydberg series of states as the ionization 
threshold is approached. 



To help our analysis, we list qualitatively different 
sources of error in the predictions from any TDDFT cal- 
culation. We refer to them as the four deadly sins: 

• The sin of the ground-state: Errors in the un- 
derlying ground-state DFT calculation. If the KS 
orbital energies are wrong to begin with, TDDFT 
corrections cannot produce accurate results. 

• The sin of locality: Errors due to local (or 
gradient-corrected) approximations to an adiabatic 
/ xc (rr'), i.e., properties that require non-locality in 
|r-r'|. 

• The sin of forgetfulness: Phenomena missing 
when the adiabatic approximation is made, i.e., 
properties that require non-locality in time, i.e., 
memory. 

• The sin of the wavefunction: Even if the exact 
u xc (ri) is used, solution of the TDKS equations 
only yields the TDKS non-interacting wavefunc- 
tion. The true wavefunction may differ so dramat- 
ically from the KS wavefunction that observables 
evaluated on the latter may be inaccurate. 

The remainder of this essay is a discussion of the vari- 
ous areas of TDDFT applications and development, and 
the challenges presently facing us. We begin in the mid- 
dle, with the linear response regime, where most of the 
applications presently are, then go to non-perturbative 
phenomena, and end with ground-state applications. 

We start with applications to the excitations of 
atoms and molecules. An important point we wish 
to emphasize here is the wealth of prediction made by 
any TDDFT approximation. The simplest real system of 
interacting electrons is the He atom, and even it has a 



rich and complex optical absorption spectrum. Return- 
ing to Fig. we note that a calculation of the optical 
spectrum of the bare ground-state KS system, often [lsj 
looks quite similar to the exact one, with TDDFT cor- 
rections merely shifting and resizing peaks. In Fig. 
we zoom out a little, and sec the ionization threshold at 
lu = I and the infinite Rydberg series of excitations just 
to its left. If one calculates the optical response of N 
non-interacting electrons in the exact KS ground-state 
potential, i.e., what we call the KS response, its ioniza- 
tion threshold is in exactly the right place, by virtue of 
the DFT version of Koopman's theorem|74j|. 

From the very earliest calculations of transition 
frea uencies |l 61 1 1 7| . it was recognized that the inaccuracy 
of standard density functional approximations (LDA, 
GGA, hybrids) for the ground-state XC potential leads 
to inaccurate KS eigenvalues. Because the approximate 
KS potentials have incorrect asymptotic behavior 
(they decay exponentially, instead of as — 1/r), the KS 
orbital eigenvalues are insufficiently negative, the ioniza- 
tion threshold is far too low, and Rydberg states are often 
unbound. This is therefore a ground-state sin. 

Given this disastrous behavior, many methods have 
been developed to asymptotically correct potentials [Til 
Frjj ]. Any corrections to the ground-state potential are 
dissatisfying, however, as the resulting potential is not a 
functional derivative of an energy functional. Even mix- 
ing one approximation for v xc (r) and another for / xc has 
become popular. A more satisfying route is to use the 
optimized effective potential (OEP) method [7^ Fif and 
include exact exchange or other self-interaction-free func- 
tionals. This produces a far more accurate KS potential, 
with the correct asymptotic behavior. The chief error is 
simply the correlation contribution to the position of the 
HOMO, i.e, a small shift. All the main features below 
and just above I are retained. 
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Why has the poor quality of ground-state potentials 
not impeded the rapid growth of TDDFT calculations for 
excitations in quantum chemistry? For many molecules, 
the lowest excitations are not Rydberg in character, and 
the orbitals do not depend on the large-r behavior of the 
potential. But there are important cases where the prob- 
lem does show up. The 'fruitfly' of TDDFT benchmarks 
is the 7r — > 7r* transition in benzene. This occurs at about 
5 eV in a ground-state LDA calculation, and ALDA shifts 
it correctly to about 7 eV[l3- Unfortunately, this is in 
the LDA continuum, which starts at about 6.5 eV! So 
how is it possible to get this right in ALDA? 

The answer is that ALDA usually yields good oscillator 
strengths, even for states pushed into the continuum [ifij]. 
The reason is simple, and was suggested long ago in e arly 
photoabsorption calculations by Zangwill and Soven|8l|. 
The KS LDA potential looks very much like the exact one 
(especially in the interior, as the occupied orbitals yield a 
good approximation to the true density), shifted up by a 
constant, due to the lack of derivative discontinuity |7J|. 
The shift pushes the Rydberg states into the continuum, 
but retains their contribution to the optical spectrum. 
Likewise for the benzene transition. Hence ALDA can 
still be used and trusted for that transition. Moral: Just 
because it's in the continuum, doesn't mean it's not right. 

The cartoon of Fig. |5] changes when bonds are 
stretched. Of particular interest in biochemistry 
are charge-transfer excitations, especially between 
weakly bonded molecules. Capturing these seems un- 
natural within TDDFT, for the simple reason that the 
numerator in x m Eq. J5J vanishes as the molecular 
separation goes to infinity. Thus, in the density-density 
response, the oscillator strength for these transitions is 
exponentially small. Recently 82], it has been shown how 
to build an empirical approximation to an adiabatic / xc 
that can capture these effects, but it is one that grows 
exponentially with |r — r'|. Thus this is a sin of locality, 
in which local approximations to / xc miss a qualitative 
feature. 

A presently o pen qu estion is the extraction of dou- 
ble excitations 83, 84] . In the adiabatic approximation, 
these are lost from the linear response equations. When 
a double excitation lies close to a single excitation, ele- 
mentary quantum mechanics shows that / xc must have 
a strong frequency dependence |if|- Thus this problem 
is due to the adiabatic approximation: a sin of forgetful- 
ness. Recently, a post-adiabatic TDDFT methodology 
has been developed[iE EH f° r including a double exci- 
tation when it is close to an optically-active single exci- 
tation, and works well for small dienes[82|]. It had been 
hoped that, by going beyond linear response, non-trivial 
double excitations would be naturally included in, e.g., 
TDLDA, but it has recently been proven that, in the 
higher-order response in TDLDA, the double excitations 
occur simply at the sum of single-excitations. Thus we do 
not currently know how best to approximate these exci- 



tations. This problem is particularly severe for quantum 
dots, where the external potential is parabolic, leading to 
multiple near degeneracies between levels of excitation. 
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FIG. 4: Same as Fig. El but now including the second ioniza- 
tion threshold 



Concerns about both ionization potentials and double 
excitations are combined when we consider more of the 
optical response of the He atom. Zooming out just a little 
more in Fig. 01 we see that there is of course a second 
ionization potential, when a second electron is stripped 
off the atom or molecule. For reference, in a He atom, 
I = 24.2 eV, and I 2 = 54 eV. But the bare KS response 
contains only the first threshold. It has no structure at 
all in the region of the second ionization. Our simple 
density functional approximations to / xc tend to shift 
and resize peak positions. It is very difficult to imagine 
density or orbital functional approximations to / xc that 
can build in the w-dependence needed to create a second 
threshold. 

Lastly in this section, we mention recent progress 
in developing a theory for electron scattering from 
molecules. This was one of the original motivations for 
developing TDDFT. One approach would be to evolve a 
wavepacket using the TDKS equations, but a more di- 
rect approach has been developed[88j, in terms of the 
response function \ of the N + 1 electron system (as- 
suming it is bound). This uses similar technology to the 
discrete transition case. Initial results for the simplest 
case, electron scattering from He + , suggest a level of ac- 
curacy comparable to bound-bound transitions, at least 
for low energies (the most difficult case for traditional 
methods, due to bound-free correlation|89]). 

A key question that often arises is the need for time- 
dependent current DFT (TDCDFT), or not. The 
Runge-Gross theorem proceeds by first proving a one-to- 
one correspondence between currents and scalar poten- 
tials. Obviously, the current is needed in the presence of 
time-dependent magnetic fields, but in their absence, is 
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it necessary? By continuity, dn/dt = — V • j(rf), so that 
the density is uniquely determined by the current (up 
to its initial value), but not vice versa. It would seem 
preferable to stay within the simpler density functional 
theory where possible. A careful examination of the con- 
ditions of applicability of the Runge-Gross theorem to 
finite systems |90j shows that all atoms and molecules sat- 
isfy the necessary conditions of potentials vanishing suf- 
ficiently rapidly as r — ► oo. However, early work showed 
that the gradient expansion (the origin of GGA's for the 
ground state) fails within TDDFT, but behaves reason- 
ably within the current theory, yielding the Vignale-Kohn 
approximation [oil l92j ] for the response kernel, which has 
frequency dependence. 

These questions become relevant to the optical re- 
sponse of bulk insulators. The Dyson- like Eq. Q 
becomes a matrix equation with indices of the recipro- 
cal lattice vectors G for each perturbation of wavevector 
q. As q — > 0, to find the optical response, any local 
approximation to / xc produces a negligible correction 
to the RPA response (f xc = 0), as the Hartree contri- 
bution (correctly) blows up as 1/q 2 - Thus, to have a 
noticable effect, the XC kernel must have a 1/q 2 compo- 
nent as q — > 0. While this effect is sometimes referred 
to as 'ultra'-non-local, we prefer to call it simply non- 
local, as the range of non-locality is precisely that of the 
Hartree contribution. The optical response of the solid 
can be found within TDDFT by perturbing the system 
with a long- wavelength perturbation of wavevector q, and 
by carefully taking q — > 0. This requires extending the 
RG theorem to periodic Hamiltonians 93] . On the other 
hand, a q — calculation, with just the period of the 
lattice, is possible within TDCDFT, and the non-local 
contribution in TDDFT appears as a local contribution 
within TDCDFT, with no unusual non-locality needed 
in the current density. For example, the VK approxi- 
mation produces a finite correction, whereas LDA and 
GGA do not. Thus, in this and other cases, TDCDFT is 
not strictly necessary, but provides a more direct descrip- 
tion and a route to extract quantities that are non-local 
in TDDFT. Similar observations apply to the polariz- 
abilities of long organic polymers. ALDA and GGA 
greatly overestimate these quantities, but VK often does 
much better. 

In this context, an important challenge is the proper 
description of excitonic peaks in the o ptic al s pect ra of 
insulators. It was recently demonstrated [94 l95l l96 | that 
with complicated orbital-dependent approximations for 
/ xc , which were derived from the Bethe-Salpeter equa- 
tion, excitonic effects can be described perfectly. How- 
ever, the presently available schemes require a GW cal- 
culation in the first place. There remains the challenge to 
find sufficiently simple (possibly current-dependent) ap- 
proximations that are able to describe excitonic effects. 

Another interesting question in the optical response of 
insulators is the one of the gap. It is well-known from 



ground-state DFT that the gap in the spectrum of KS 
eigenvalues (the KS gap) differs from the true gap by a 
quantity called the derivative discontinuity 3] . Ignoring 
excitons within the gap, shouldn't TDDFT correct the 
KS gap to yield the true gap? The answer is, yes, but 
again the XC kernel that does this must be very sophis- 
ticated, just as in our double ionization example. Since 
Xs develops an imaginary part for frequencies above the 
KS gap, the kernel must have a branch cut that exactly 
suppresses this in order to widen the gap. There is a 
close analogy to the problem of charge-transfer excita- 
tions: Remove an electron from the donor to infinity. 
This costs the ionization energy /don- Then move the 
electron from infinity to the acceptor. In this way one 
gains the energy — Aacc- So the excitation energy is 

Ac 1 j A _ c LUMO c HOMO , A ( Acc ) „4,™. Q 

iXa — iDON - ^ACC — £aCC — DON + ' wnerc 



(ACC) 



is the discontinuity in the ground-state w X c( r ) of 
the acceptor. This formula is reminiscent of the band 
gap in insulators 3J. Furthermore, the next-order correc- 
tion is the Coulomb-interaction between the electron on 
the acceptor and the hole on the donor which, in solids, 
corresponds to the exciton binding energy. 

On the boundary between extended systems and 
molecules is transport through single molecules con- 
nected to bulk metal leads |97|. There is enormous inter- 
est in this as a key component in future nanotechnology. 
Present formulations use ground-state density function- 
al to describe the stationary non-equilibrium current- 
carrying state But several recent suggestions con- 
side r this as a time-dependent problem |64l 199 . llOOL llOll 
102], and use TD(C)DFT for a full description of the sit- 
uation. Only time will tell if TDDFT is really needed 
for an accurate description of these devices. In the spe- 
cial case of weak bias, XC corrections to the Landauer 
formula are missed by local approximations, the sin of 
locality |l03j . 

Next we turn our attention to beyond-perturbative 
regimes. Due to advances in laser technology over the 
past decade, many experiments are now possible in 
regimes where the laser field is stronger than the nu- 
clear attraction. There are a whole host of phenomena 
that TDDFT might be able to predict: high harmonic 
generation, multi-photon ionization, above-threshold ion- 
ization, above-threshold dissociation, etc. For high 
harmonic generation, TDDFT c alcul ations have been 
rath er s uccesful for atoms .104, 105] and molecules 
[lil llOfij . In the near future, this might become very 
important for the generation of atto-second laser pulses 
For multi-photon ionization, the rela- 



tive proportion of double to single ionization for He, while 
given much better in approximate TDDFT calculations 
than in previous calculations assuming a sequential mech- 
anism, still does not show the sam e pronoun ced features 
(the 'knee') seen in experiments |ll0t 111 l| . The elec- 
tron spectra from above-threshold ionization have 
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recently been calculated within TDDFT [H[TT3. Since 
the ionization yields and above-threshold ionization spec- 
tra depend on probabilities extracted directly from the 
wavefunction, these errors are suspected to be sins of the 
wavefunction, mentioned above. An important task for 
the future will be the design of more realistic expressions 
for ionization probabilities or, more generally, transition 
probabilities as functionals of the time-dependent den- 
sity or the time-dependent KS or bitals . First steps in 
this direction can be found in Ref. jllflj - 

While the need for more accurate KS potentials was 
first noticeable in calculating excitations, it is even more 
acute in strong laser fields. To ensure an atom or 
molecule ionizes accurately in an approximate TDKS cal- 
culation, Koopmans' theorem (I = —€homo) should be 
well-satisfied, and this again requires using OEP exact 
exchange jzj or other self- interaction- free function- 
als. 

The field of quantum control has, so far, mainly con- 
centrated on manipulating the motion of the nuclear wave 
packet on a given set of precalculated potential energy 
surfaces, the ultimate go al be ing the femto-second con- 
trol of chemical reactions |l!3| . With atto-second pulses 
available, the control of electronic dynamics has come 
within reach. A marriage of optimal-control theory with 
TDDFT appears to be the ideal theoretical tool to tackle 
this situation. However, it will bring with it its own diffi- 
culties and challenges for approximate functionals. Using 
the functi onal algorithms developed by Rabitz and others 
|l!4illl5| . we can find the optimal pulse that drives a He 
atom from its ground state to its first excited state, ls2p. 
(In practice, we do not reach exactly 100% occupation, 
due to a finite penalty factor) . 

Now repeat this experiment on non-interacting elec- 
trons sitting in the same potential. Such a pulse cannot 
be found, i.e., the non- interacting system is not control- 
lable, whereas the interacting system is. The two non- 
interacting electrons must follow the same time evolution 
as they start from the same initial Is orbital and are ex- 
posed to the same laser field. Hence the time-dependent 
wave function of the two non-interacting electrons must 
have the form: 

$(ri,o-i,r 2 ,<72,t) = ^(ri,t)v?(r 2 ,t)xs(CTi,cr 2 ) (11) 

where xsi&i,^) represents the (antisymmetric) spin- 
singlet part of the wavefunction. But we want to maxi- 
mize the occupation 



|< $(T)|$ ls , 2p >| 



(12) 



of the time-propagated wavefunction $(T) at the end, T, 
of the laser pulse in the lowest excited state 

$is,2 P = -^=(<^i s (ri)^ 2p (r 2 ) + ^ 2 p(ri)^i 5 .(r 2 ))x5(CTi,cr 2 ) 

(13) 



Expanding the final wavefuntion <j>(T) in the complete set 
of single-particle orbitals representing the eigenfunctions 
of the unperturbed system, one sees that the best possible 
occupation is 50%. 

If TDDFT is used to describe the multi-electron dy- 
namics, how would one properly define the control target, 
i.e., the functional to be maximized? Choosing, as con- 
trol target, the overlap with an excited-state Kohn-Sham 
determinant does not seem to be a good idea in view of 
the above dilemma. If, on the other hand, the time- 
dependent density of a fully controlled (Is) 2 to (ls,2p) 
transition of the interacting system is given, would an ex- 
act TDKS calculation reproduce such an optimal evolu- 
tion? The answer is yes, but v xc (rt) m ust be very special 
to do so. To see this, take the density evolution from the 
exact Schrodinger equation, and invert the TDKS equa- 
tion for the single (doubly-occupied) time-dependent or- 
bital. The final state KS potential is very odd, producing 
the density of two orbitals of different symmetry from a 
single doubly-occupied orbital [llfij . 

For a subset of cases in which molecules are exposed 
to strong fields, the nuclear motion can be treated clas- 
sically. The electrons then feel the Coulomb field of clas- 
sically moving nuclei as well as the laser field. In this 
case, the electronic motion is well described by ordinary 
TDDFT. However, when nuclear trajectories split, 
e.g., when a molecule has a 50% chance of dissociation 
in a given laser pulse, t he classical tre atment fails. A 
multicomponent TDDFT [111 Ell has been devel- 
oped for electrons and nuclei which, in principle, handles 
such a situation. In practice, one needs to develop appro- 
priate approximations that can build in the appropriate 
physics of, e.g., dissociating nuclei. Again, in this ca se, 
orbital-dependent functionals appear crucial |ll7illl9| |. 

Finally, and fondly, we return to the ground state. The 
last general application mentioned was the odd-sounding 
process of using TDDFT to generate ground-state 
approximations, via Eq. By inserting an approx- 
imation for / xc , we get an approximation to E xc . Most 
importantly for biological systems, Eq. provides a 
route to van der Waals forces for separated pieces of 
matter, and so is being much studied by developers. In 
particular, the coefficient in the decay of the energy be- 
tween two such pieces (Cg in E — > —Cq/R 6 , where R is 
their separation) can be accurately (within about 20%) 
evaluated using a local approximation to the frequency- 
dependent oolarizabilitv [2ll fl2d Il2ll Il22j . Recent work 
shows that the response functions of TDDFT can yield 
extremely accurate dispersion energies of monomers |l23| . 
More recently, the frequency integral in Eq. © has been 
done approximately, yielding an explicit non-local den- 
sity functionally applicable at all separations. 

One can go t he other w ay, and try using Eq. for 
all bond lengths ^24 Il25| . Such calculations are costly, 
as the functional is very high up on Jacob's ladder of 
sophistication, including both occupied and unoccupied 
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KS orbitals 63]. However, they have the merit of being 
entirely non-empirical and, where successful, can be used 
as a starting point for new approximations. In fact, Eq 
@ provides a KS density functional that allows bond- 
breaking without artifical symmetry breaking [22|. In 
the paradigm case of the H2 molecule, the binding energy 
curve has no Coulson-Fischer point, and the dissociation 
occurs correctly to two isolated H atoms. Unfortunately, 
simple approximations, while yielding correct results near 
equilibrium and at infinity, produce an unphysical repul- 
sion at large but finite separations. This can be traced 
back [2^ to the lack of double excitations in any adiabatic 

/xc- 

We end with a summary. Most importantly, TDDFT 
has become extremely popular as a method for calculat- 
ing electronic excited-state energies in chemistry. In this 
arena, it has become as robust (or as flaky, depending 
on your perspective) as ground-state DFT, and is being 
used to really understand the photochemistry of many 
interesting problems. Furthermore, our favorite crude 
density functional approximations from the ground-state 
serve well here. We are very thankful for this, and it 
has led to tremendous interest in further methodological 
development. 

In principle, however, TDDFT yields predictions for 
an enormous variety of phenomena, and electronic excita- 
tions are only the tip of an iceberg. We have mentioned a 
few. Even limiting ourselves to linear response, there are 
double excitations, second ionization thresholds, optical 
response of solids, gaps in solids, transport through single 
molecules. Combining with the fluctuation-dissipation 
theorem, TDDFT yields a route to van der Waals forces 
and bond-breaking with symmetry problems. In strong 
fields, there are high harmonic generation, multi-photon 
ionization, above-threshold ionization, quantum control, 
and quantum nuclear motion. 

For some of these areas, simple application of density 
functionals within the adiabatic approximation, works 
well, but for many, such methods miss some qualitative 
features (e.g., double excitations, or non-locality in the 
response of solids). A now standard step upward in so- 
phistication is to use orbital-dependent functionals (at 
least, among developers), and these cure some of the dif- 
ficulties (e.g., the first ionization threshold or the polar- 
izability of long-chain polymers). But such functionals 
are unlikely to cure all of the problems (e.g., inclusion 
of double excitations or defining the target in quantum 
control) for properties that are of interest experimentally 
and technologically. We happily look forward to many 
interesting years of development to come. 
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